Current induced spin flip scattering at interfaces in noncollinear magnetic multilayers 
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We show that when one drives a charge current across noncollinear magnetic layers the two 
electron Coulomb scattering creates a spin flip potential at interfaces. This scattering is found when 
the interface potential is updated due to the spin accumulation attendant to charge flow, and it 
contributes in linear response to the current. With this scattering there is an injection of transverse 
spin distributions in a layer that propagate, and that in the steady state lead to spin currents 
transverse to the magnetization. 
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Our current knowledge of the electronic structure of 
magnetic multilayers is based on effective single electron 
approximations. In nonmagnetic (normal) layers states 
are spin degenerate while in the magnetic they are spin 
split; therefore one can construct a coherent combina- 
tion of spin up and down states in the normal layers, 
while in the magnetic this is not possible because there 
is no coherence (correlations) between states in spin split 
bands i This loss of spin coherence as one crosses a nor- 
mal/ferromagnetic layer (N/F) interface has been identi- 
fied as the origin for the inability of injecting transverse 
spin currents from a normal metal into a 3d transition- 
metal ferromagnet. 2 The studies to date of the coherence 
transmitted by interface scattering have concluded that 
spin currents across this interface can induce only a trans- 
verse coherence between electrons on different sheets of 
the spin split Fermi surface of the ferromagnet. 2 Here 
we revisit this question and find a current driven cor- 
rection to the transmission amplitudes coming from the 
two electron Coulomb interaction at interfaces that al- 
lows the transverse spin coherence to be carried by the 
electron distribution function for each spin split band. As 
we show while this is an out of equilibrium effect, its in- 
clusion in transport calculations enters in linear response. 
In another publication we show that this current induced 
coherence make it possible for spin currents to be contin- 
uous across the interface in the steady stated Inasmuch 
as one includes current driven accumulations in the bulk 
of the layers when calculating the steady state currents 
we conclude that to discuss spin currents in noncollinear 
magnetic structures it is necessary to self consistently 
update the transmission and reflection coefficients due 
to the presence of out of equilibrium spin distributions 
at interfaces. 

To relate distribution functions f ss '(k,r) across an in- 
terface one resorts to the transfer matrbsi 

f>,(k,0 + )= J2 T mm , 8S ,f> m ,(k',e k ,0-) 

mm'k' 

+ £ R mm , ss ,f< m ,(k",e k ,0 + ), (1) 

mm'k" 

where the transmission and reflection coefficients (prob- 



abilities) are 

Tmm'^ss' — T^mm' ' ss' — s * ^m' s' i 

Rmm' =?ss' — R-mm' ss' — ?™m s * ^m's'i (^) 

and t mg (k', k, e k ), r ms (k" , k, E k ) are the transmission and 
reflection amplitudes for elastic scattering. At the inter- 
face between a normal and ferromagnetic layer the mo- 
menta k,k',k" are either ku or k m , for a majority and 
minority bands of the ferromagnetic layer and k n for a 
normal metal layer. We make the "isotropy" assumption 
that allows one to write one coefficient for all momen- 
tum directions, albeit one should still keep a band index, 
i.e., T mm i^, ss '(k',k,e k ) ~ T mm >=> S3 >(e k )£ In a spin po- 
larized single electron description the spin split bands of 
the magnetic layer have only one spin index, s = s', for 
each momentum k ,i.e., f ss ' = fs, however when one is 
propagating currents across noncollinear structures one 
has to envisage the possibility of off diagonal components 
of the distribution function f^ s ,(k, + ) with s ^ s' within 
a screening length about the interface. This represents 
a coherence between states in different bands; it does 
not represent an admixture of spin within the band al- 
though the scattering we find is capable of this. D The 
indices to, to' refer to the spin of the electron in the un- 
split Fermi sea of the normal metal layer; in nonmagnetic 
layers it is possible to construct coherent combinations 
of up and down spin states to ^ to'. Scattering at inter- 
faces in metallic structures is localized within a screening 
length about the interface and represented by amplitudes 
t mt ->m r ,i"mi -^m[ ] the first describe the transmission from, 
say, the left layer to the right, the second the reflection 
back into the left layer, and the m's are the components 
of the electron's spin. 

The scattering potential, t matrix, for conduction elec- 
trons at an interface between two metals has contribu- 
tions from differences in, the periodic background po- 
tentials, the kinetic energy and the Coulomb interaction 
between the electrons; when written as an effective one 
electron potential only the latter changes due to current 
driven spin accumulation on the Fermi surface. To find 
these corrections we focus on the electron-electron con- 
tribution to the t matrix. The one electron scattering 
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amplitudes at interfaces are derived form the two elec- 
tron Coulomb interaction 

Vcoulomb = 2 X! V ( k isk2s'k3 s 'k4 s )cl iS cl 2S ,Ck 3 s'Ck i s, 

S,8 

where V{k\ s ..k4 :S ) are the matrix elements of the 
Coulomb interaction between states on either side of the 
interface. For our purposes the range of integration is 
limited to the screening length about the interface. As 
we are working at the interface between dissimilar well 
screened metals the translational invariance of the back- 
ground is broken; therefore we do not reduce the four 
to three momenta. To determine the one electron scat- 
tering amplitude, f mi _ >mr ,r mi _ >m ', that arises from this 
Coulomb interaction we write the Coulomb interaction 
between electrons in the one electron states in the bulk 
of the layers on either side of the interface, and to reduce 
it to a one electron operator we take the expectation 
value over a pair of annihilation and creation operators 
in Eq.@. In the lowest order distorted wave Born ap- 
proximation the spin dependent part of the correction to 
the transition matrix between one electron states on the 
two sides of the interface is& 

t m ^s = (k 2 ,s \t op \ k n ,m) , (4) 

where 

top = ~7;g{£k 2 )^2 ( C !im C fe3s) y(klmk2sk3 s k nm )cl zS C knm , 
kik 3 

(5) 

\k n ,m > refers to states on the Fermi surface of the nor- 
mal layer, and < fca, s| to those on the ferromagnetic spin 
split Fermi surface, so that k 2 — ku and k m . For elastic 
scattering 5(ek n — £k 2 ) the density of states g{£k 2 ) enters 
when one averages the amplitude t m ^ s over the energy 
£k n - While the states in the ferromagnetic layers are 
pure spin states, e.g., only an up spin goes in the major- 
ity sheet, we are looking for current induced coherences 
in the distribution function, a statistical density matrix, 
between states of opposite spin; therefore we do not im- 
mediately associate a spin s to a state k 2 , e.g., we will be 
looking for matrix elements for both spin directions on 
each sheet of the Fermi surface of the ferromagnetic layer 
that are induced when the system is out of equilibrium. 

The expectation value (c\ im Ck 3 s) is between states 

on opposite sides of the interface and is found by self 
consistently evaluating it there; as we will be interested 
only in elastic scattering the energy of the states entering 
this expectation value are the same, £k t — £k 3 , but their 
directions in k space can be different. At the interface 
between normal and ferromagnetic layers (N/F) the scat- 
tering potential, while spin dependent, is diagonal in spin 
space, i.e., there is a unique spin direction, and there are 
no elements m ^ s when the system is in equilibrium, and 

4 ira Cfc 3S ) ~ n l s nt S sm ; it also follows that = t s S sm . 



When driving current across a magnetic multilayer spin 
accumulates on the Fermi surface of the normal layers 
so as to support a spin current across them in steady 
state, therefore in the presence of a spin current the spin 
polarized distribution, 5n s) of the Fermi surface of the 
magnetic layer M; upstream from the N/F r interface is 
superimposed on the Fermi surface of the normal layer. 

For noncollinear magnetic layers the magnitude of 
the current driven accumulation in the interfacial region 
N/F r coming from the left magnetic layer, 8n" Lt , is un- 
certain; this can only be ascertained by a calculation of 
the transmission amplitudes at N/F r when one superim- 
poses 8n l s on the Fermi surface (at the Fermi energy) of 
the normal layer i In linear response we do know it is 
quantized along the magnetization Mi of the magnetic 
layer upstream while the operators entering the expecta- 
tion value in the interfacial region ( see Eq.JSJ) are quan- 
tized along M r the magnetization of the ferromagnet at 
the N/F r interface. To determine the additional contri- 
bution from 6n l s to the transmission amplitude, Eq.JjjJ, 
it is necessary to rotate the accumulation referred to Mi 
to states quantized along M r ; we find 

&top = A ma (k2 1 k n )c\ 2S Ck n m, (6) 

where 

A m s(k 2 , k n ) = V s (k s k n ; £k 2 ) 
f [Sn\ nt cos 2 9/2 + Sn mt s sin 2 9/2] S sm \ 
\ -| [Snf * - 6n mt s ] sin 6S s> _ m J ' [ ' 

V s (k 27 k n ;£ k2 ) = -^g s {£ k2 ) ((V(k 2s ,k nm ;£ F )}} , (8) 

((V s (k 2s , k nm ;s))) = 
j^pp J dfl kl J dQ, k3 V s (k 1 ,k 2s ,k 3 ,k nm ), (9) 

Here <5rt™* is the occupancy of the spin states on the 
Fermi surface at the interface which have been polarized 
by the upstream magnetic layer Mi. 

The out of equilibrium accumulation and current 
transmitted across an interface N/F by the equilib- 
rium T matrix is found from Eq. JJ| by integrating 
over all k states. For the diagonal spin components 
fss' = fT ^ss' the out of equilibrium distribution only 
exists on the Fermi surface, because /^C 1 '^!:,^) ~ 
<5 fm{k')5 mm i5{£k — £f)- For spin transport across non- 
collinear magnetic layers the state \k n ,m > is quantized 
along the magnetization Mi of the magnetic layer up- 
stream while the operators in the t matrix are quantized 
along M r the magnetization of the ferromagnet at the 
N/F interface, i.e., 9 = cos~ 1 ( Mi ■ M r ), therefore it is 
necessary write | k n , m > in terms of \k 2l s >. By ro- 
tating this state we find the average of the transmission 
coefficient in equilibrium is 

T mm ^ss{£F) = N 2 { cos2 9/25 

sm "h sin 6 /IS s—m } , 

(10) 
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where t s = t^ s S sm . As the spin diagonal term ~ S sm in 
Eq.Q is proportional to the accumulation or current it 
does not enter in linear response in the T matrix that 
connects the spin diagonal part of out of equilibrium dis- 
tribution functions, f> and 8f m (k'), as they themselves 
produce the accumulation and current when integrated 
over the Fermi surface. However the off diagonal (spin 
flip) term~ <5 s _ m in Eq.Q is uncompensated because the 
accumulation coming from the right layer at the N/F r is 
diagonal in spin space, and cannot offset an off diagonal 
term coming from the left. At the other interface of the 
normal layer, N/F/ , the accumulation arising from the 
right layer 8n r s cannot be fully compensated from that 
arising from the left layer and we have an uncompen- 
sated contribution 5t op ~ <5 s ._ m . 8 

The current driven coherences f^ 8 , with s ^ s' do 
not represent populations, that in equilibrium are given 
by the Fermi Dirac distribution function, therefore their 
contribution to the transmission across an interface is 
not limited to the Fermi surface (sk ^ £f) and has con- 
tributions from the entire Fermi sea. As the scattering 
potential, Eq.JfJJ, is proportional to the accumulation or 
current its contribution in linear response to the spin 
coherence transmitted across an interface is found by 
using the equilibrium distribution in the normal layer 
/ (£fc,0 - )£ mTO / in Eq.JTJ, and by using an equilibrium 
transmission amplitude to find the update to the equilib- 
rium T matrix in the presence of a current. As the dis- 
tribution function entering the right hand side of Eq. 1(1)1 
does not depend on the spin index m, we find, in linear 
response, the contribution of the spin flip transmission 
amplitude Eq.© to the off diagonal transmission coeffi- 
cient ST mm ^, ss r for states on the Fermi surface comes 
from integrating over the Fermi sea 

/>/(Mf,0+) 

- E fe , ro {C * K ls - + Mm* * C*< } 0-) 

= 6n? t anO{$£Vf[<rvr] aa > +&* [a Xr } ss ,} , (11) 
where 



t s V* = J de k f (e k )gn(£k) ((t s (k2,k n ;Sk)V*{k s k n ;e k ))) , 

. (12) 

£s(fe, k n ; Sk) is the equilibrium transmission amplitude, 
and the angular brackets are defined in Eq.©. Here the 
out of equilibrium densities are Sn z = i [6rif — 5n>i\, and 
to arrive at this result we made the isotropy assumption 
for the transmission coefficients so that the angular av- 
erages over these coefficients are done independently of 
those over the distribution function. 

This transmission coefficient is current driven so that 
in linear response it acts on the equilibrium distribution 
function at the Fermi surface in Eq. © to produce an out 
of equilibrium distribution f^", ; whereas the transmis- 
sion Eq. ( HOfl acts on the out of equilibrium distribution 



on the Fermi surface. By writing Sn" lt = a5n l z and 

S < = \ E [^Jmm / dn k ,5f m (k% (13) 
m 

we find the two contributions, Eas. l|lU|) and i|ll|) . to the 
T matrix can be written as 

f>,(k,£ F ,Q + ) = ^2T mm ^ ss (e F )S ss ,Sf m (k') + iasin6> 

m m 
X [<J Zl ] mm {&* [<Ty r ] ss , + &tJ?[<Ts r ] aa ,}Sf m &). 

(14) 

Whereas the first term is the equilibrium scattering at 
the interface in the conventional approach, the second 
depends on the spin accumulation or current at the inter- 
face which is proportional to 5f m {k'), and produces an off 
diagonal component of the spinor density matrix at the 
interface. While V s is a sum over the Fermi surface ( see 
Eq.© with e — sf) and t s over the entire Fermi sea their 
difference is made up in the integration over the equilib- 
rium distribution function that enters Eci. (|12|) which is 
over Fermi sea. We can only do a meaningful compar- 
ison, of the transmission amplitudes in equilibrium and 
another when one superimposes 5n z on the Fermi surface 
of the normal layer, by calculating the interface scatter- 
ing for these two situations; interalia this will determine 
the constant a which indicates how much of the spin ac- 
cumulation 6n l z penetrates into the interfacial region so 
as to affect the interface scattering potential^ 

By solving the Boltzmann equations of motion we have 
found that the magnitude of the "out of equilibrium" ST 
relative to T controls the amount of transverse spin accu- 
mulation, but its very existence guarantees the continu- 
ity of spin currents at the interface when a steady state 
is achieved, i.e., to eventually achieve a continuous spin 
current in the steady state it is necessary that the inter- 
face scattering potential has matrix elements that inject 
transverse spin distributions into a magnetic layer with 
well defined momentum so that they can propagate past 
the interface, and into the bulk of a layer. 3 Without this 
additional scattering at the interface we have found a dis- 
continuity in the spin current at all times; while this is 
found by setting the explicit time derivative of the distri- 
bution function, dtf(k,r,t), to zero this solution differs 
from steady state inasmuch as there is a constant loss of 
transverse spin current at the interface, i.e., it is ques- 
tionable whether one can call this a true steady state. 
The transverse spin current j y and j' x (this exists if the 
transmission amplitudes are complex) we calculate in the 
ferromagnetic layer, by using Ea. (|llfl or 1)14(1. arises from 
injecting the transverse component of the incoming spin 
current from a normal layer in such a manner that it 
excites a transverse mode of propagation in the bulk of 
the ferromagnetic layer. Our scattering potential pro- 
vides for the transfer of the transverse spin distribution 
across an interface in such a manner that it can prop- 
agate past the interfacial region and into the bulk of a 
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magnetic layer; it do not suffer from the dephasing pre- 
viously found when one considered the transverse spin 
current that arises from scattering onto different sheets 
of the Fermi surface, e.g., 

T m m^.ss> ~ [t s {kM)t*s'{k m )] {i/2 sin 0<5 S /_ S } , (15) 

where fc^z/m represent the Fermi momenta of the ma- 
jority/minority bands. Among other things the distri- 
bution function f> s , (k, ep, + ) found by using this T 
in Eq.Q does not have a unique velocity v(k),as ku ^ 
k m , therefore it cannot be used to calculate a current 
within the Boltzmann approach. The existence of this 
scattering between bands limits transverse currents to a 
region 1-2 monolayers of the interface^ but it does not 
negate the current induced interface scattering mecha- 
nisms, Ea. (|ll(l . that allow for injection into a coherent 
mode of transverse spin propagation on another length 
scale. 

The spin current transmitted by the conventional coef- 
ficients, Eq. frroj). is parallel to the local magnetization; 
it does not have a transverse component. If the interface 
scattering Eq. Q is unable to scatter electrons into trans- 
verse spin distribution functions f ss > (k) s' ^ s which can 
propagate past the interfacial region when the current 
is first turned on, the transverse component of the spin 
current never leaves the interfacial region; not even in 
steady state. In a manner of speaking the transverse 
spin accumulation will remain confined to the interface 
and one can talk about the spin current as being dis- 
continuous even in the steady state. The two electron 
scattering across an interface Eq.JJJ) creates the scatter- 
ing potential, Eqs.@ and (|llfl that allows for a coherent 
transmission of spin information at the interface. How- 
ever this only exists at the interface, within a screening 
length of the interface; without a current this electron 
distribution does not go beyond the interface scattering 



region. In the presence of a current once created at the 
interface the transverse distribution propagates accord- 
ing to the Boltzmann equation; as the distribution f ss i 
for s' ^ s does not commute with the Hamiltonian in 
a magnetic layer ( which is described in the spin polar- 
ized single electron approximation) we show in another 
publication this leads to distributions that precess, and 
that it is the exchange splitting in this approximation 
J(k) that controls this precession. 3 Eventuallythis leads 
to spin accumulation transverse to the magnetization of 
the magnetic layer, so that in steady state the spin cur- 
rent is continuous across the N/F interfaced 

In conclusion our scattering potential provides for the 
transfer of the transverse spin distribution across an in- 
terface in such a manner that it can propagate past the 
interfacial region and into the bulk of a magnetic layer. 
The spin flip scattering potential Eq.© lies outside the 
conventional spin polarized single electron treatment of 
the scattering at a N/F interface as one does not usually 
envisage when calculating the transmission amplitude a 
distribution from the spin polarized Fermi surface of a 
neighboring magnetic layer M; being superimposed onto 
the normal layer £ However when one uses an approxi- 
mate conductivity that is local (short ranged), e.g., when 
using only the bubble conductivity in the Kubo formal- 
ism or the Boltzmann layer by layer approach, in order 
to account for the long range nature of the conductivity 
we are required to posit this in the presence of an electric 
fieldiiSii The strength of the spin flip potential relative 
to the ordinary scattering at the interface can be ascer- 
tained from calculations of the transmission amplitudes 
in the presence of spin accumulation in the normal layer. 
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